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Context. We consider a simple self-gravitating disk, made of two fluid components characterized by different effective thermal 
speeds and interacting with one another only through gravity; two-component models of this type have often been considered 
in order to estimate the impact of the cold interstellar medium on gravitational instabilities in star-dominated galaxy disks. 
' Aims. This simple model allows us to produce a unified description of instabilities in non-viscous self-gravitating disks, some 
f****) , originating from Jeans collapse, and others from the relative motion between the two components. In particular, the model 
suggests that the small streaming velocity between the two components associated with the so-called asymmetric drift may be 
the origin of instability for suitable non-axisymmetric perturbations. 

Methods. The result is obtained by examining the properties of a local, linear dispersion relation for tightly wound density waves 
© ' in such two-component model. The parameters characterizing the equilibrium model and the related dispersion relation allow 
us to recover as natural limits the cases, known in the literature, in which the relative drift between the two components is 
<yH • ignored. 

' Results. Dynamically, the instability is similar to (although gentler than) that known to affect counter-rotating disks. However, 
^» , in contrast to the instability induced by counter-rotation, which is a relatively rare phenomenon, the mechanism discussed in 
V~j ■ this paper is likely to be rather common in nature. 

Conclusions. We briefly indicate some consequences of the instability on the evolution of galaxy disks and possible applications 
^ , to other astrophysical systems, in particular to protostellar disks and accretion disks. 
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1. Introduction 

As clearly noted in investigations of the problem of spiral 
structure in galaxies, the dynamics of galaxy disks de- 
pends significantly on the mutual interaction between gas 
and stars. For the general issues related to the properties 
of waves and instabilities, such gas-star interactions have 
often been studied by referring to a very simple, idealized 
two-component model in which one distinguishes a "hot" 
component (meant to represent mostly the older stellar 
disk) from a "cold" component (representative of atomic 
and molecular hydrogen and young stars) interacting with 
each other only through gravity. In practice, calculations 
have been carried out by treating one component by means 
of the equations of stellar dynamics and the other by fluid 
equations, or, for simplicity, by treating both components 
by means of simple barotropic fluid models characterized 
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by different effective thermal speeds; other investigations 
have implemented these concepts by means of suitable N- 
body simulations. It has thus been ascertained that even 
small amounts of gas can have significant influence on the 
Jeans stability of galaxy disks (e.g., see Lin & Shu 1966 
Lvnden-Bell ITTM7I Graham ITMTl Miller, Prendergast & 
QuirklP3701 QuirkHnm Kato lT$731 Jog & Solomon lP3Mal 
I1984bl Sellwood & Carlberg lT$gl Bertin & Romeo 113551 
hereafter BR88). 

Similar models have also been used to examine situa- 
tions in which the two components are in relative motion 
with respect to one another. Among the earliest investiga- 
tions, we may recall those by Sweet |T9j53J ). Lynden-Bcll 
(1967), Wilson & Lynden-Bell (unpublished; as quoted 
in Lvnden-Bell ITMTjl . and Marochnik & Suchkov ifTMH 
I1974JI . Often in such early studies the "Population II" 
stars are treated as a non-rotating component, in a ki- 
netic model, and, for simplicity, the geometry adopted 
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is that of a homogeneous three-dimensional system (see 
also the extensive discussion of the interaction between 
gas and stars in the monograph by Saslaw I1985fl . Later 
on, similar studies were further motivated by the issue of 
disk- halo interaction (Mark 1 976f) : so, in the latter inves- 
tigations one component represents the rotating disk and 
the other the non-rotating or slowly rotating halo mate- 
rial. These analyses then brought out the possibility of a 
kind of two-stream or Kelvin-Hclmholtz instability, thus 
beyond Jeans instability, with applications to the ampli- 
fication of bending waves in relation to the problem of 
galaxy warps (Bertin & Mark 1980S Bertin & Casertano 
IT5£2l revisited by Nelson & Tremaine ITO) . With sim- 
ilar tools, others have explored the possibility of inter- 
actions between different gaseous components in galaxies 
(see Nelson E7H| Waxman HWfM I1979bfl . 

More recently, renewed attention has been drawn to 
the study of instabilities excited in disks in relative mo- 
tion (Sellwood & Merritt 1994; Lovelace, Jore & Haynes 
119971 hereafter LJH97; Hunter, Whitaker & Lovelace 
Comins et al. IT337I Mandica l2TOT|l by the discov- 
ery of the curious phenomenon of counter-rotation, noted 
in a number of disk galaxies (Rubin, Graham & Kenncy 
AMU Rix et al.Em Merrifield k Kuiiken lT§9H Bertola 
et al.QjSEI Jore, Broeils & Havnes HMHI Prada et al.QjUl 
Thakar et aL UMTI Bettoni et al. lM)T|l . 

Counter-rotation appears to be a relatively rare, and 
possibly short-lived phenomenon. Thus, such relative mo- 
tion between two disk components is likely to be only 
rarely involved in the dynamics of galaxies. In contrast, 
in this paper we focus our attention on a potential source 
of instability that might be ubiquitous: because of the so- 
called asymmetric drift associated with the finite velocity 
dispersion of the hot component (e.g., see Bottema ll99,S(l . 
any disk of stars and gas would in general be character- 
ized by small, but finite relative motion between the two 
components. 

In this paper we will show that, in principle, the small 
velocity difference between two components associated 
with the asymmetric drift is a source of instabilities in 
the disk. This result will be obtained by inspection of a 
local, linear dispersion relation for tightly wound density 
waves derived from a simple two-fluid model in Sect. 2. In 
Sect. 3 we will study the conditions for the onset of the 
asymmetric-drift instability and its main characteristics. 
Finally, in Sect. 4 we will describe possible applications to 
different astrophysical contexts and draw the main conclu- 
sions of our paper. Given the fact that such asymmetric- 
drift velocity is small, for applications to galaxy disks, this 
paper is meant to offer an investigation preliminary to a 
fluid-kinetic analysis that is planned for the future. 

While this work was being completed, we were re- 
minded of an interesting study made in the context of 
the dynamics of protostellar disks ( Weidenschilling 1 1 977|l . 
about the viscous drag exerted by the gas of the protostel- 
lar disk on the population of planetesimals, which presents 
significant analogies with the mechanism addressed in this 



paper. We will briefly comment on this point in the con- 
clusion section. 

2. The fluid model and the relevant dispersion 
relation 

In a rotating self-gravitating axisymmetric fluid disk at 
equilibrium the radial gravitational force is balanced by 
rotation, with a contribution from the pressure gradient: 



where f2, cr, and p are the equilibrium angular velocity, 
surface density, and pressure of the thin disk, and <E> is the 
total gravitational potential. For cool disks, the pressure 
gradient is small and generally neglected. Therefore, many 
studies of the dynamics of such disks simply refer to the 
angular frequency as defined by the relation rft^ = d^/dr 
and to the related epicyclic frequency defined by Kq = 
4fig[l + {l/2)dhxfl /d\nr}. 

We now consider an infinitesimally thin, self- 
gravitating, two-component disk. In general, the two com- 
ponents are characterized by different equivalent acoustic 
speeds Cj and surface densities <7£. We label the two com- 
ponents by the subscripts c and h (in relation to their 
equivalent acoustic speeds, so that c c < c/j). Similarly 
to other previous stability analyses, we assume that the 
only interaction between the components occurs via the 
gravitational field. In contrast with previous analyses that 
focused on the Jeans instability and assigned the same 
equilibrium angular velocity f2o to the two components 
(e.g., see BR88), in this paper we allow for the presence of 
relative motion between the two components, so that in 
general f2 c ^ Qh, and, in turn, for a difference in the asso- 
ciated epicyclic frequencies k c ^ Kh- The special limiting 
case in which c c — > Ch and f2 c — ► — Qh thus corresponds 
to the case of counter-rotating components, a situation 
already discussed in the literature, as mentioned in the 
Introduction. We argue that, because of the different size 
of the relevant asymmetric drift dictated by Eq. the 
angular velocities of the two components would in general 
be (although slightly) different even when the two compo- 
nents are not counter-rotating. 

2.1. Applicability of the two-fluid model 

The two interpenetrated non- viscous fluids studied in this 
paper are to be taken as the simplest model in which 
we can include several key factors (temperature, density, 
and velocity differences, within the non-trivial geometry 
of an inhomogenous thin layer) that determine the stabil- 
ity of a two-component self-gravitating disk. Some impor- 
tant features, such as finite thickness, present in real self- 
gravitating disks are not included. Therefore, for any spe- 
cific application (for example, to the dynamics of galaxy 
disks) it will be necessary to check to what extent the ba- 
sic approximations implicit in this model are justified and, 
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if necessary, to extend the calculations to a more realistic 
(and thus more complicated) model. On the other hand, 
studying this simple model has the advantage of leading 
to a unified picture of a variety of instabilities that may 
affect a self-gravitating disk. 

One particular point of concern that will emerge from 
the following analysis is related to the applicability of the 
two-fluid model to cases in which the relative speed be- 
tween the two components is small. In fact, it is well known 
that, for a kinetic system made of two components in rel- 
ative streaming, the fluid description becomes inadequate 
when the relative speed is small (e.g., see Lynden-Bell 
IPMTl Krall & TrievelpieceEHZIi- The instability noted in a 
two-fluid description may remain, but is usually modified 
because the underlying mechanism of instability is asso- 
ciated with an inverse Landau damping (i.e., the effects 
of wave-particle resonances) that only a kinetic descrip- 
tion can properly describe. Thus, in view of applications 
to a disk made of stars and gas (a system that is nei- 
ther fully kinetic nor fully fluid), the present study should 
be completed by investigating the mechanisms in a mixed 
fluid-kinetic model; this extension will be addressed in a 
separate paper. 

2.2. Local dispersion relation 

We now consider tightly wound linear density perturba- 
tions of the form an = &u exp[i(— u)t+m6+ J kdr)] under 
the WKB ordering m/(r\k\) = 0(e), with the epicyclic 
expansion parameter defined as e = c; l /(rK/ l ). We will 
take the azimuthal number m to be positive, so that trail- 
ing disturbances are characterized by radial wave number 
k > 0. To lowest order, the density response for each com- 
ponent (resulting from the relevant continuity and Euler 
equations) is given by 

k 2 $i 

Each component experiences the joint perturbed gravita- 
tional potential <I>i, which is calculated from the Poisson 
equation 

- 1*1$! = 27tG((7 1c + a lh ) . (3) 

By combining the above equations to eliminate $i (for a 
complete derivation see Cava l2U(1Hl . we obtain the disper- 
sion relation, which we write in dimensionless form as 

v A - 2i W 3 - A^v 2 + 2rjA^v + A (0) = . (4) 

Here the relevant coefficients AW = (\k\; Q\, a, /?, A) 
for n ^ 1 and A« = A^(\k\;Q 2 h ) are defined as 

A® = 1 + A^ + A - a\k\ + ^Ql/3\k\ 2 , (5) 
A^=l-\k\ + ^Ql\k\ 2 , (6) 



A® = (1 + A)A« - a\k\ + -^Ql(3\k\ 2 - (7) 
l Q l(a + /3)\k\ 3 + ^PQ 4 h \k\ 4 • 

Here k — 2nGahk/ 'k 2 and v = (ui — mfi/j)//^ are 
the dimensionless radial wavenumber and the dimension- 
less Doppler-shifted frequency, respectively. The quantity 
Qh = ChKh/(jrGah) is the axisymmetric stability parame- 
ter relative to the hot component. 

Thus, with respect to the one-component stability 
analysis (which is uniquely controlled by the axisymmet- 
ric stability parameter), the dispersion relation depends 
on four additional dimensionless parameters a, /3, A, and 
rj. For the first two (positive definite) parameters we follow 
the definition of BR88, i.e. a = cr c /(Jh and (3 = c 2 /c\ < 1 
are the relative density and "temperature" ratios. The 
two new parameters are defined in terms of the angu- 
lar velocity ratio S = fl c /flh (which, for simplicity, was 
taken to be constant in the derivation of Eq. ifljl) as 
rj = m(S—l)(Clh/Kh) and A = 5 2 — r) 2 — l; therefore, in the 
limit £l c — > we have A — ► and 77 — > 0. Note that 
r\ vanishes for axisymmetric perturbations (m = 0) even 
when the two components are in relative motion (S ^ 1). 
Since we ascribe the difference in angular velocity to the 
different weight of the pressure term in Eq. l(T|l. we expect 
the cooler component to rotate faster, so that in general 
|I|>1. 

2.3. Recovery of known limiting cases 

For (5=1 the two components are corotating exactly, 
which is the case considered by BR88, while for 6 = — 1 
the components are counter-rotating, which corresponds 
to the case considered by LJH97 (who, at variance with 
the present two-fluid description, considered a mixed fluid- 
kinetic model and a kinetic-kinetic model). 

The one-component dispersion relation v 2 = 
A^> (\k\; Qf^ is recovered by taking the limit (3 — > 1, 
6 — > 1 (so that -q — ► and A — ► 0), and by letting 
a become vanishingly small; the second, odd branch (see 

Sect. 2.2 in BR88) of the dispersion relation then reduces 

- 2 

to the condition v 2 = 1 + (l/4)Q 2 i |fc| for modified sound 
waves in the "gas tracer". For a two-component system 
the limit of exactly corotating fluids is obtained by letting 
rj = and A = so that the dispersion relation reduces 
to V A ~A {2) V 2 +A { ^ =0, with4 0) (l*l;Qfc,«,/3) = 
a|fc| + (l/4)Q2 /3 |fc| 2 _ ( i / 4 ) g2 (a + / j ) |fc| 3 + (1/16)/ 3Q4|^ 4 . 

thus the marginal stability curve at the basis of the anal- 
ysis of BR88 (see their Eq. (14)) is recovered by setting 

4 0) = 0. 



G. Bertin and A. Cava: Local stability of self-gravitating fluid disks 



4 

3. The asymmetric-drift instability 

3.1. Axisym metric disturbances: reduction to standard 
Jeans stability conditions 

We first consider the case of axisymmetric disturbances, 
that is the case of r\ = 0, with A = A = S 2 — 1 > 0. Here 
the dispersion relation reduces to v 4 — A^v 2 + = 0. 
Therefore, the marginal stability condition is given by 
= 0. From Eq. J7J) it is clear that, even if A ^ 0, 
the condition is the same as for the A = case, provided 
the relevant parameters a and are suitably rescaled. In 
fact, if we define a 1 = a/(l + Ao) and 0' = 0/(1 + Ao), we 
find that the condition A^ (\k\;Ql, a, 0, A ) = is equiv- 
alent to the condition A^ (\k\; Q\, a', 0') — 0. Therefore, 
for axisymmetric disturbances the stability condition is 
the same as for the standard Jeans instability in two- 
component disks. In particular, local stability at all wave- 
lengths is guaranteed when Qh > Q(a',0'), which re- 
places the local axisymmetric stability condition for one- 
component fluid disks Q > 1; the function Q = Q(a,0) 
is illustrated in Fig. [5] of BR88. Note that, in general, for 
A > we have Q(a\0') < Q(a,0), which corresponds 
to a stabilizing effect on m = perturbations. We should 
also point out that the rescaling of a and (3 that has al- 
lowed us to make use of the marginal stability condition 
described in BR88 does not correspond to exactly map- 
ping the roots of the dispersion relation with 5 > 1 into 
those of the case with 5 = 1. 

3.2. Instability for non-axisymmetric perturbations at 
small values of rj 

We now refer to the case when the relative motion be- 
tween the two components is small, that is 77 <§C 1, with 
A = Ao — i] 2 < Ao, and consider an equilibrium case that 
is marginally stable with respect to the Jeans instability, 
so that Q h = Q(a/(l + A ),/?/(l + A )). The "most dan- 
gerous perturbations" are those with wavenumber close to 
I k I , defined by the relation A 0) ( | k | ; Q 2 , a/ ( 1 + A ) , 0/ ( 1 + 
Ao)) = 0; close to the quantity A*- ' is expected to 
be negative and small, because A < Ao- In addition, we 
expect A^ > (the last relation follows from the fact 
that the presence of a second component is known to have 
a destabilizing effect on axisymmetric perturbations, i.e. 
Q(a',0') > 1). As to the value of the frequency for such 
conditions close to marginal stability, we expect it to be 
also small, and thus we can approximate the relevant dis- 
persion relation as a quadratic 

-A^v 2 + 2nA^v + A^ ~0 . (8) 

Therefore, in the wavenumber range close to |fc| in which 

V 2 (A^) 2 + A^A^ <0 (9) 

we find a pair of complex-conjugate solutions one of which 
corresponds to instability, with dimensionless growth rate 



vi = 0(r)). The real part of the frequency for these com- 
plex solutions is given by 

v * ~ < V , (10) 

which, as anticipated, corresponds to a pattern speed 
Q p = ujji/m in the range Qh < Q p < 51 c . 

These facts are illustrated in Fig.^, where we compare 
the dispersion relation in the vicinity of marginal stability 
from the exact (Eq. 4) and from the "reduced" (Eq. 8) 
expression, showing that the approximate quadratic rela- 
tion can be accurate even for relatively large values of rj. 
In Fig.fJlwe also illustrate the behavior of the coefficients 
A^> and A^\ in line with the arguments given above to 
derive the onset of the instability analytically. 

3.3. Non-axisymmetric perturbations studied in the 
general case 

At this point we can move on to study the general case, 
when the new parameters A (or S) and 77 that appear 
in the dispersion relation cannot be taken to be small. 
As described previously, we have in mind two different 
kinds of situations of astrophysical interest. The case of 
counter- rotation is associated with S « — 1 so that we ex- 
pect 77 and A to have finite values. As noted earlier in the 
Introduction, this case has already been addressed in the 
literature. The second case, where the velocity difference 
of the two components, being generated by the asymmet- 
ric drift, is small, has been dealt with in the previous sub- 
section; here, in order to consider 77 = 0(1), we should 
take m to be large. In this new regime of open waves the 
basic dispersion relation Eq. would be modified and 
there would exist additional destabilizing mechanisms (see 
Bertin et al. 119891 note the threshold value indicated by 
their Eq. (5.8), which, for a flat rotation curve, reduces to 
Q 2 = 3). For this, a consistent asymptotic analysis might 
be developed in an approach similar to the one followed 
in this paper, but for simplicity we omit such analysis in 
this paper. 

In Figs . 13141 we illustrate the qualitative behavior of the 
dispersion relation at fixed value of S for varying 77 and for 
varying m (with fixed flh/^h)', in each plot the dashed and 
solid lines represent the imaginary and real part of the di- 
mensionless frequency v 1 respectively. Since, at variance 
with the case of axisymmetric perturbations, the unstable 
root has a non-vanishing real part of the frequency, we 
are in the presence of a convective instability (see LJH97). 
Following the general arguments used by LJH97, we can 
also confirm for this general case of 77 = O(l) that in- 
stability can develop when vr(vr — rj) < 0, i.e. when 

Qh ^ Qp ^ Qc 

In conclusion, the dispersion relation indicates the ex- 
istence of unstable roots with dimensionless growth-rate 
(imaginary part of v) up to 0.1 — 0.2 for small values of 
m, and formally up to 0.4 for higher values of m. 
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Fig. 1. The dispersion relation in the vicinity of marginal sta- 
bility from the exact (Eq.^J top frames) and from the "re- 
duced" (Eq.|S] bottom frames) expression, for two cases char- 
acterized by S — 1.2, rj = 0.1 and S — 1.2, rj = 0.3. In each 
frame the real (solid lines) and imaginary (dashed lines) part 
of v are given as a function of the dimensionless wavenumber 
k. The other parameters have been set to the following values: 
a = 0.1, (3 = 0.1, Q h = Q(a/P,f3/S 2 ) « 1.14. The "reduced" 
dispersion relation turns out to be accurate even for r\ — 0.3. 
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Fig. 2. Illustration of the generation of the unstable roots. 
Superposed to the dispersion relation shown in the second 
frame of Fig.Q (for the case 8 — 1.2, r\ = 0.3, with a = 0.1, 
f3 = 0.1, Qh ~ 1-14) we plot the coefficients ,4 (0) (dotted curve) 
and AP^ (dash-dotted curve) for values of k close to k. This 
diagram is given in support of the arguments that lead to the 
condition (Eq.|5J for the onset of the instability. 



3.4. Relevance to astrophysical applications 

As we will discuss further in the following Section, the 
main context for astrophysical applications that we have 
in mind is that of normal spiral galaxies. 

The strategy that we have used has been to focus on 
a condition of marginal stability with respect to axisym- 
metric perturbations in the two-component disk (Qh = 
Q(a,/3)) in order to demonstrate that the relative mo- 
tion of the two components, independently of axisymmct- 
ric Jeans-type instabilities, is destabilizing. The instabil- 
ity that we have addressed can obviously occur under 
more general conditions than Qh = Q(a 1 f3), as shown 
in Fig.|31 in particular, instability can occur for a — 0.1 
and (3 = 2.5 10 -2 (parameter values suggested by studies 
of the solar neighborhood) even for Qh > 1.5. 

The asymmetric drift instability, being driven by the 
relative motion of the two components, is expected to be 
ubiquitous and to coexist (as mentioned in the first para- 
graph of §3.3) with other sources of instability. As a re- 
sult, the destabilizing effect that we have discussed, even 
if present in realistic simulations of two-component disks, 
may be difficult to single out from other sources of insta- 
bility and from a number of effects related to the specific 
modeling used in each simulation (e.g., see Orlova et al. 
|2TjO"21 see also Noh et al. lTMTl Korchagin & Theis HMUjl . 

Interestingly, the instability described in this paper is 
expected to operate non-linearly to couple the two com- 
ponents (that is to keep cold and hot components under 
physical conditions not too far apart), counteracting the 
effects of the cooling processes present in the cold gas, 
which favor the decoupling of the two components. This 
makes the asymmetric drift instability participate in the 
" thermostat" of self-regulation (necessary for grand design 
spiral structure to be maintained; see Bertin & Romeo 
1988] Bertin et al. I1989|l . The natural emergence of self- 
regulation has been demonstrated recently in numerical 
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experiments of the dynamics of self-gravitating protostel- 
lar disks (Lodato & Rice OUMjl . 

Another point that we wish to emphasize is that the 
study that we have performed is not limited to "small- 
scale instabilities" (the complex unstable root shown in 
Figs.0]and[3]occurs at k of order unity, which does not cor- 
respond to "small scale" wavelengths); in fact, this prop- 
erty distinguishes the mechanism studied in this paper 
from the destabilizing effect on axisymmetric stability as- 
sociated with the presence of cold gas, which indeed tends 
to act on the small scale. 

4. Discussion and conclusions 

In this paper we have shown that a self-gravitating thin 
disk made of a hotter and a cooler component can be un- 
stable (with respect to non-axisymmetric density waves) 
because of the different amounts of asymmetric drift as- 
sociated with the two components, even when the disk is 
stable in the sense of Jeans (with respect to axisymmet- 
ric density waves). The unified analysis performed here is 
based on a local dispersion relation from which we can eas- 
ily recover the limits of previous stability investigations. 
In particular, we can recognize the two-component Jeans 
stability conditions, as studied by BR88, and the insta- 
bility of counter-rotating disks, as discussed by LJH97 
(in a technically more complicated study). For the case 
of small relative velocity between the two components, 
which is the focus of the present paper, and of low values 
of the azimuthal wavenumber m, the growth rate of the 
instability has been derived analytically. Here we aimed at 
investigating the asymmetric drift instability mechanism 
by itself and were not interested in assessing the stability 
or instability of a disk under specific physical conditions. 
For this latter purpose, it is well known that a number 
of other factors (such as the role of thickness; e.g., see 
Romeo IT992[I should be properly taken into account. The 
collective effects realized in specific contexts will depend 
on the conditions of the basic state and on the non-linear 
phenomena that take place during evolution. 

For applications to galaxy disks, the fluid- fluid analysis 
presented here is preliminary to a more realistic analysis 
in which the cold component is treated as a fluid, while the 
hotter component is treated in terms of the collisionless 
Boltzmann equation. Resonance effects, to be described 
by the equations of stellar dynamics, are expected to oc- 
cur when the velocity difference between the two compo- 
nents is small, i.e. when it becomes comparable to the 
velocity dispersion of the stars associated with the hotter 
component; note that this condition refers to the velocity 
dispersion in the azimuthal direction, which is smaller by 
approximately a factor of k/(20) than the radial velocity 
dispersion, which is responsible for the asymmetric drift. 
A proper description of resonances will lead to technical 
complications well beyond those addressed in the fluid- 
kinetic study of LJH97. Because of its physical origin, the 
asymmetric drift velocity is a priori bound to be relatively 
small, even for those galaxies (e.g., NGC 5064, NGC 3200, 



and NGC 2815) in which such velocity difference has been 
found to be, in magnitude, rather large (in the range of 
50 to 100 km/s; see Vega Beltran et al. I2()01l and Pizzella 
et al. l27JUl|) . 

For galaxy disks, the instability studied in this paper is 
likely to induce turbulence and heating in the cold compo- 
nent and to contribute to the generation of a net inflow of 
material toward the center. In principle, one might think 
of a possible connection with the topic of extra-planar 
gas that is currently the subject of several investigations; 
this is the discovery that many galaxy disks, in addition 
to their "normal" thin HI layer, possess an "anomalous" 
thick HI component characterized by slower rotation (e.g., 
for NGC 2403, see Fraternali et al. E002 for UGC 7321, 
see Matthews & Wood 1271(131 for NGC 4559, see Barbieri 
et al. l2005|l . However, for these systems it is likely that the 
two HI components do not interpenetrate, but rather oc- 
cupy different heights with respect to the equatorial plane; 
in this case, if instability occurs, the mechanism might be 
different from the one studied in this paper, and more in 
line with processes of the type investigated by Waxman 
l|l97MlTi79b|l . 

Finding a mechanism able to act as a source of turbu- 
lence in disks brings us to consider possible applications 
to the long-standing issue of finding an adequate expla- 
nation for the anomalous viscosity that is thought to op- 
erate in accretion disks (see Shakura & Sunyaev ©73 ) . 
The present mechanism would contribute in those cases 
in which the various physical processes support the coex- 
istence of two components at different levels of rotation 
around the central object. As a kind of two-stream insta- 
bility, it might operate even for plasma disks. 

This latter application naturally brings us to consider 
the case of protostellar disks. In fact, the mechanism stud- 
ied in this paper appears to be the collective counter- 
part to the "viscous drag mechanism" (Weidenschilling 
11977(1 that initially appeared to act against the picture 
of planet growth from planetesimals. Recent simulations 
have shown that such drag may actually have a beneficial 
role in the direction of planet formation (see Rice et al. 
2004|. This paper suggests that a specific collective pro- 
cess, independent of viscosity, may operate and contribute 
in a way similar to that of mechanisms present in viscous 
disks (see also Youdin & Goodman 2005 ) . The very simple 
analytical study provided in the present paper offers a de- 
scription, at the linear level, of the collective instabilities 
underlying such interesting process. 
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Fig. 3. The dispersion relation for varying rj. In each frame 
the real (solid lines) and imaginary (dashed lines) part of v 
are given as a function of the dimensionless wavenumber k, for 
77 = 0,77 = 1,77 = 2,77 = 3 (from top to bottom). The other 
parameters have been set to the following values: <5 = 1.2, 
a = 0.1, (3 = 0.1, Q h = Q(a/S 2 ,(3/8 2 ) « 1.14. 
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Fig. 4. The dispersion relation for varying m. In each frame 
the real (solid lines) and imaginary (dashed lines) part of v 
are given as a function of the dimensionless wavenumber k, 
for m = 0, m = 1, m = 2, m = 3 (from top to bottom). 
The other parameters have been set to the following values: 
5 = 1.2, a = 0.1, (3 = 0.1, Q h = Q(a/S 2 ,(3/8 2 ) w 1.14, with 
the additional condition of flat rotation curve Q. h /K h = \/2/2. 
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Fig. 5. The asymmetric drift instability for a disk stable with 
respect to axisymmetric instabilities, Qh > Q(a,@). Here a — 
Q.1,/3 — 2.5 10 -2 , values suggested by conditions in the disk 
of the Milky Way. The vertical axes in the second and fourth 
frames start from -0.50 to better illustrate the imaginary part 
of the root of the dispersion relation (Eq.^J. 



